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Abstract 

1^ , We prove the following theorems: 1) The Laurent expansions in e of the Gauss 



hypergeometric functions 2Fi{Ii + ae, h + be; I3 + ^ + ce; z), 2-^1 (A + ^ + ae,l2 + 



2 + te; /a + 2 + ce; z) and 2Fi{Ii + ^+ae,l2 + be; /s + f + ce; z), where Ii, /2, h,P, Q 
^ . are arbitrary integers, a, b, c are arbitrary numbers and e is an infinitesimal pa- 

r^ I rameter, are expressible in terms of multiple polylogarithms of g-roots of unity 

with coefficients that are ratios of polynomials; 2) The Laurent expansion of the 



Q ■ Gauss hypergeometric function 2Fi{Ii + ^ + ae, I2 + be;!^ + ce; z) is expressible 



in terms of multiple polylogarithms of q'-roots of unity times powers of logarithm 



^^ I with coefficients that are ratios of polynomials; 3) The multiple inverse rational 

—». ' sums YlTLi—? — — \^Sai{j — 1) ■ ■ ■ Sa {j — 1) and the multiple rational sums 

> 

X 

H , - ^ - ■ 

5t ! 4) The generalized hypergeometric functions pFp^i{A + ae; B + be,^ + Bp-i; z) and 

Fp^i{A + ae,-+Ap;B + be;z) are expressible in terms of multiple polylogarithms 



Ejli r(i+j) F*^"!^-^' "!)■■■ 'S'ap(j - 1)> where Sa{j) = J2i=i p: is a harmonic series 
and c is an arbitrary integer, are expressible in terms of multiple polylogarithms; 



p-t p 



Q 



with coefficients that are ratios of polynomials. 
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1 Introduction: Feynman diagrams and special func- 
tions 

High-precision theoretical predictions for physics at the LHC and the ILC demand the 
inclusion of higher-order radiative corrections. The results of perturbative calculation are 
expressible in terms of Feynman integrals [JLj. However, in order to obtain physical results, 
it is necessary to construct the Laurent expansions of Feynman diagrams about the integer 
value of the space-time dimension [2] (typically d = A — 2e). For the parametrization of the 
coefficients of such e expansions, a lot of new functions have been introduced by physicists 
during the last few years [3]IHl5]|6] . Some of these new functions are also generated in a 
different branch of mathematics [7f8|9llT0] . At present, it is unclear if there is some 
limitation on the types of functions generated by Feynman diagrams or if the "zoo" of 
new functions is an artifact of the techniques used. In particular, the statement that the 
results of such calculations can be written in terms of a restricted set of special functions 
will allow one to use a restricted set of programs for the numerical evaluation of physical 
results. Another application is related to the evaluation of so-called single-scale diagrams, 
where an explicit prediction of possible transcendental constants can be done [TT] . 

The strategy of such a kind of analysis is well know in the theory of special functions 
and the analytical theory of differential equations [12] . As is well known, any Feynman 
diagram satisfies a system of linear differential or difference equations with polynomial 
coefficients [T3lll4lll5fl6lll7j . In modern mathematical language, such a system can be 
associated with the Gelfand-Karpanov-Zelevinskii functions or Z)-modules [18]. So, any 
question regarding the zoo of special functions generated by the e expansion of Feyn- 
man diagrams could be reduced to the problem of constructing Laurent expansions of 
D-modules (hypergeometric functions [19]) about certain values of their parameters. Un- 
fortunately, a unique hypergeometric representation of Feynman diagrams besides the 
so-called a representation [1] does not exist. Using the latter representation, it has been 
shown recently that, for single-scale diagrams, i.e. diagrams where all kinematic variables 
are proportional to each other so that one of them can be factored out, all coefficients 
of the e expansion can be understood, up to some normalisation factor, as periods in 
the Kontsevich-Zagier formulation ^20j. Another useful representation, which is closely 
related to the corresponding property of generalized hypergeometric functions, is the 
Mellin-Barnes representation of Feynman diagrams [21]. Using the Mellin-Barnes repre- 
sentation, it is possible to write the result in each order of e in terms of multiple sums, 
which sometimes can be expressed in terms of special functions [22] • Since the power of 
a propagator is integer in covariant gauge and any (irreducible) numerator is expressible 
in terms of an integral of the same topology with a shifted power, which is again inte- 
ger [T^IITS] . it is enough to only consider hypergeometric functions of several variables 
with integer values of parameters. (In general, the number of variables is equal to the 
number of kinematic invariants minus one.) Fortunately, when some of the kinematic 
invariants are proportional (or equal) to each other, the number of variables in the proper 
hypergeometric series can be reduced. But the price of this reduction is the appearance 
of rational values of parameters. All known exactly solvable cases [23ll24ll25f26f27j have 



confirmed this observation. Typically, only integer and half-integer values of parameters 
are generated, and only recently inverse cubic values have been discovered |27j . 

Recently, there has been essential progress in understanding what type of functions 
are generated by the e expansion of hypergeometric functions. Besides the pioneering con- 
struction of the e expansions of hypergeometric functions [28] using harmonic series [29] 
or so-called multiple (inverse) binomial sums [30|31|32j . there are now a few independent 
techniques for the construction of analytical coefficients in the e expansions of hypergeo- 
metric functions about integer and half-integer values of parameters and the summing of 
multiple series [32|33f34f35f36ll37f38f39f40f41j . However, the extension of these results to 
the case of rational values of parameters is still a mystery. There is just one publication 
[36] devoted to the analysis of e expansions of hypergeometric functions about a special 
configuration of rational parameters, the so-called, "zero-balance" case. Specifically, two 
types of sums have been analyzed in Ref. [36], namely in Eqs. (51) and (62). To our knowl- 
edge, they correspond to the sums covered by our Theorems A and B below, respectively, 
in the case c > 1. 

The aim of the present paper is to derive an algorithm for the construction of all-order 
e expansions of generalized hypergeometric functions pFp_i and multiple (inverse) rational 
sums. The present consideration is based on appropriate extensions of the generating- 
function approach [32|42|l39j and the differential-equation technique [37|38f40j to the case 
of rational values of parameters. 

In particular, we will prove the following theorems: 

• Theorem I: 

The all- order e expansions of the Gauss hypergeometric functions 
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Fi (Ji + as, h + fe; /a + 2 + ce; z) , (la) 
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2F1 (/i + f + ae, h + he; h + ce]z) , (lb) 

2F1 (/i + f + ae, h + 6£; Is + I + ce; z) , (Ic) 

2F1 (/i + f + as, /2 + H + 65; /3 + £ + ce- z) , (Id) 

where {Ik} are integers, a, b, c are arbitrary numbers, and e is an arbitrarily small 
parameter, are expressible in terms of multiple polylogarithms (or multiple polylog- 
arithms times powers of logarithm) with coefficients that are ratios of polynomials 
with complex coefficients. 

Theorem A 

The multiple inverse rational sums 



E 



r(j)r (1 



r(i+i-f)^ 



s<.,(i-i)---s.,(i-i), (2) 



where Sa{j) is a harmonic series, defined as Sa{j) = J2k=i I^' ^'^^ ^ '^^ '^''^V iiT'teger, 
are expressible in terms of multiple polylogarithms with (i) complex coefficients if c > 



1; and (ii) with coefficients that are ratios of polynomials with complex coefficients 
ifc<0. 

• Theorem B 

The multiple rational sums 



prr(i + j)r(i + f)^^ 

where Sa{j) is a harmonic series, defined as Sa{j) = J2k=i 'h' ^^^ ^ ^'^ '^''^y ^'^^c^'er, 
are expressible in terms of multiple polylogarithms times powers of logarithm with (i) 
complex coefficients if c> 1; and (ii) with coefficients that are ratios of polynomials 
with complex coefficients if c < 0. 

Theorem C 

The all-order e expansion of the generalized hypergeometric functions 

pFp„i (^A+ae, s + Ja; B + he- z^ , (4a) 

Fp_i (A+as] B + he, f+Ji; z\ , (4b) 



p-t p-i , -- . "-, - ■ "-' g 



where A, B are lists of integers and Ii, I2 are integers, are expressible in terms of 
multiple polylogarithms (or multiple polylogarithms times powers of logarithms) with 
coefficients that are ratios of polynomials with complex coefficients. 

This paper is organized as follows. In Section [2], we will prove Theorem I. In Sec- 
tion [31 we will present an analysis of multiple (inverse) rational sums and prove Theorem 
A and Theorem B. In Section HJ the results of Theorem A and Theorem B will be ap- 
plied to hypergeometric functions to prove Theorem C. In Section [5l we will demonstrate 
that, for physically interesting kinematics, the two-loop sunset diagrams are expressible in 
terms of generalized hypergeometric functions with quarter values of parameters. Appen- 
dices \M and [B] contain some basic information about multiple polylogarithms, which are 
a particular class of hyperlogarithms, and the iterative solution of systems of differential 
equations. 

2 Gauss hypergeometric function: notations 

The Gauss hypergeometric function [13] uj{z) = 2Fi{a,b;c; z) can be defined as the so- 
lution of a second-order differential equation of Fuchsian class [12J with three regular 
singular points at 2; = 0, 1, cxd, as 

z— + b\uj{z), (5) 
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f d 


\ 




( d 


— j 


z-r -hc- 


-1 


u{z) = 


2, 


dz 


l dz 


y 


V dz 



and admits the series representation about z = 0, 






fc=0 



where {a)k = T{a + k)/r{a) is the Pochhammer symboL It is the only solution analytic 
at z = and with value 1 there. 

It is well known that, between any three Gauss hypergeometric functions with the 
same argument z and parameters a, b, c differing by integers, there is an algebraic relation 
with polynomial coefficients [H], namely 



Pi(a, 6, c, z)2Fi{a + mi, b + ni;c + ki; z) + P2{a, b, c, z)2Fi{a + m2, b + n2;c + ^2; z) 
+P3(a, b, c, ^)2-Fi(a + ma, 6 + ns; c + /ca; 2;) = , (7) 

where rrij, Uj, kj G Z (j = f , 2, 3). Taking ms = n^ = k^ = and m2 = n2 = /i;2 = 1, we 
obtain, for example, by using the algorithm described in Ref. [55], 



P(a, b, c, z)2Fi{a + mi, b + Ui; c + ki, z) 



Qi{a, b, c, z)—- + Q2{a, b, c, z) 
az 



2Fi{a,b;c;z 



, ^j 1^, A^y , 



where a, b, c, are any fixed numbers and P,Qi,Q2 are polynomial in the parameters a, b, c 
and the argument z. We call the functions of r.h.s. of Eq. ([8]) basis functions and their 
first derivatives. Consequently, in order to prove Theorem A, all-order e expansions 
have to be constructed for basis hypergeometric functions. 

2.1 Differential equation approach for construction of e expan- 
sion 

2.1.1 Gauss hypergeometric functions 



Let us consider as the basis the Gauss hypergeometric function with the following set 
of parameters: u{z) 
differential equation 



of parameters: u{z) = 2Fi{— + aie, — + 02^; I — — + ce; z) . It is the solution of the 



d pi \ f d P2 \ / \ d f d p3 \ , \ /N 

z—- -\ \- aie ] z—- -\ \- a2e ] uj{z) = -— z- h ce uj{z) , (9) 

dz qi J \ dz q2 J dz \ dz q^ J 

with boundary conditions iu{0) = 1 and z^u{z)\^^q = 0. Due to the analyticity of the 
Gauss hypergeometric function with respect to its parameters, Eq. ([9]) is valid in each 
order of e, i.e. it holds for every coefficient function ujk{z) in the expansion 

CX) 

cj{z) = J2Mz)e'. (10) 

A:=0 



The boundary conditions for the coefficient functions are 

ujk{z) = (^ < 0) , (11a) 

c^fc(O) = (A; > 1) , (lib) 

^tM^)\.=o = ^^ (^^0)- (lie) 

Equation (Q can be rewritten in terms of the coefficients functions Uk as 



r. ^d fpi P2\ lp3 

dz \qi q2j zq3 



d pip2 

dz qiq2 



= (ai + a2--) z-—ujk-i + I ai — + 02— 1 uJk-i+aia2UJk-2 ■ (12) 

V zJ dz \ q2 qij 

The main idea of the approach developed in Refs. PSPU] is to find a new parametrization, 
through change of variable z — >• ^{z), and to define new functions PkiO^ related to the 
ffist derivative of the original functions uJkiO^ ^s 

so that Eq. flT2l) can be rewritten as a system of linear differential equations of ffist order 
with rational coefficients, as 

iM0=Pk{0Y.7^^ (13a) 

^PkiO = P.-i(0 E 7%- + ^^^-1(0 E J^ + ^^^-2(0 E J^^ ' (13b) 

where Aj, Bj, Cj, Dj, aj, Pj, •jj, ctj G C (j = 1, 2, ■ ■ ■ ). Then, the iterative solution of this 
system can be constructed as explained in Appendix [Bl Under the condition ujo{z) = 1 
{po{z) = 0) this solution is expressible in terms of hyperlogarithms (see Appendix |A]) de- 
pending on the parameters aj, Pj, ■jj, aj, possibly times powers of logarithm. For example, 
the ffist iteration of Eq. ( ]13bl) produces 

Pi(0 = E^^- iGi{ir,a^)) - Gi(7,;e(0))] (7, ^ 0) , (14a) 

j 

u;i(O = $^A,A{[Gi,i(a„7fc;e(^))-Gi,i(a„7fc;e(0))] (Mb) 

- [Gi(«,;e(^))Gi(7fc;e(0)) - Gi(«,; e(0))Gi(7fc; ^(0))]} (7fc,«, ^ 0) . 

The main problem is to find a general algorithm for constructing this parametrization. 
We are not able to proof that we found a solution of this problem for all possible values 
of the parameters. But for some special set of parameters, the solution is found. 

This algorithm can be applied to construct the all-order e expansion of an arbitrary 
system of differential equations of the ffist order, in particular to the generalized hyper- 
geometric function, as was done in Ref. [ID]. 

6 



2.1.2 One lower parameter is a rational number 

Let us consider the particular case that the basis function is of the form 

P 



ujiz) = 2F1 ( aiE, aiE\ 1 V ce\z 

q 



(15) 



where we assume that p, g > and p < q. The differential equation (Q takes the form 



;i-.)-f_iP 

az z q 



Z—UJk[Z) 

az 



with the boundary condition 



Let us introduce a new variable ^jj 



fli + 02 ) z--uJk-i{z) + aia2UJk-2{z) 

zJ az 



u)q{z) = 1 . 



(16) 



(17) 



z-1 



1/9 



so that 



l-^g ' 



1- z 



1 



1-z 



d 



dz q dt. 



'X-z 



\-i'i' z e 

d pi 1 fl-^'i 



dz q z q \ ^'' 

We introduce new functions pk{C,) via the differential equation 



^ir' 



d 

Z--UJk[z) 

dz 



q di 



The boundary conditions for the new coefficient functions are 

Pfc(O) =0 (A; > 0) . 



(19) 



(20) 



(21) 



Equation (TT6l) can be rewritten in terms of new coefficients functions ^^k{C) ^-^id pk{C) ^^ 
a system of two first-order differential equations, as 



1 d S,^~^ 

qdt, 1 — t,^ 



1 d 

— -77pk{0 
qdi 



(ai + 02) 



jg-l 



i-e^ i 



(22a) 
Pfc-i(0 + Q1Q2 ^^g ^fc-2(0 • (22b) 



^The proposition to use the variable ^ for the evaluation of multiple inverse rational sums was made 



in Ref. [36j . For the particular value q — 2, the variable ^ is equivalent to the variable y — — ^ 



considered in Ref. |32| . due to the invariance of the Remiddi-Vermaseren functions [3] with respect to 
the transformation z — > tx^. For q = 1 
Remiddi-Vermaseren functions in Ref. Hw 



the transformation z — > jxf • For 9 = 2, the variable S, was also applied to the parametrization of 



Now, we are in a position to proof the following result. 
Lemma I: 

The all-order e expansion of the Gauss hypergeometric function 2F1 ( cti^, 02^; l — - + c€;z 
is expressible in terms of multiple polylogarithms with arguments that are powers of q-roots 
of unity and the variable ^ defined by Eq. ( I7g|) . 
This may be written symbolically as 

2F1 ( aie,a2e] 1 Vce\z 

00 

J, s 

1 < {jm} < q 

EI=1 S^ = J 

(23) 

where s = {si, ■ ■ ■ , s^} is a multi-index and vjg are numerical coefficients {vjg G C). In 
particular, the following statement is valid: 
Corollary I: 

The analytical coefficient of e^ in the expansion of 2F1 ( ^i^, Ci2^] \ — - + ce; z\ includes 

only multiple polylogarithms of weight k with numerical coefficients. 

Proof: 

Firstly, it is necessary to show that the system of equations fl22|) can be rewritten in the 

form of Eq. fll3bp . This follows from the standard decomposition relation 



1 - ^' = n (1 - K^) ' (24) 



where we have introduced the primitive g-root of unity. 



Using the decompositioro 



Xq = exp ( i — ) . (25) 






2 



For completeness, we present a few particular cases: 

j = l AJ j=l A| j = l A^ 



where < r < g — 1, the system of equations (l22l) can be rewritten for an arbitrary value 
of p (0 < p < g — 1) in the desired form (113b|) . as 



d 



\y^ 



3=1 ^ K 



di 

d 



Pkio 



(ai + aa)^-— ^ 



1 cq 






The solution of system ( 1271) has the form 

^« dt 



^k 



7=1 -^0 



—Pk{i) 



(27a) 



(28a) 



A^o 



Pfc(0 



(ai + as)^ / - — --cq — Pfc-i(t) + aia2 V Af / p^fc-2(t) • 



The first coefficients of the e expansion are zero, 
The first nontrivial terms correspond to A; = 2Jj 



P2(0 

aia2 
aia2 



J] A^^nn (1-A^;e) = - E A^Lii (A^;e) 



ii=i 
g 



ii=i 



5^ A(^-^^Hii,i (A^,-^s A;^e) • 



ii,«i=i 



^This is equivalent to the following representation for the hypergeonietric function: 



^.^, M;2-^;z =^ 
q-p \ q J \z~l 



f fiT^)" t1-P-^ 



fi - 1 



■dt , 



(28b) 
(29) 

(30a) 
(30b) 



where z ^ \. 



Higher-order terms can be generated by iteration: 

(l\(l2 ■^""^ ■^""^ 

ii,i2=i ii=i 

i^ = {a. + a^) J2 A(^-^^)aii,i,i (A^;-^% Af-% A;^e) 

0102 . ^~^ -, 

J1J2,«1=1 

q 



+cq J2 A(^-^^)^Lii,2(A^;-^\A;^0 , (31b) 

ii,n=i 

^^ = - E ^r [(«i+«2)'Lii,i,i (A^;-^^ Af -^3^ Af + c^g'Lis (A^;0] 

01^2 • — 

-cg(ai + a2) ^ A^;4 Lii,2 (A^,^-^% Af ^ + Li2,i (A^,-^^ Af C) 

jl,J2 = l K 

q 

+aia2 E A(^■^+^■^-^^Hil,l,l(A^,-^^A^l-^^Afe) • (31c) 

iij2,n=i 

Let us apply the mathematical induction. Let us assume that Lemma I is valid up to 
order j, so that 

^.(0 = E ^/,s-Li^- iK~''^ Af -^^ ■ ■ ■ , A>— >, A>e) , (32a) 

J, s 

1 < {jm} < q 

Er 
i=l Si=J 

P.(0 = E "^Lifc (A^^ Ai-'3^ ■ ■ ■ , A^— '-, A^-e) • (32b) 

k,t 
1 < {la} < q 
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Substituting these expressions in Eq. (]28bp . we obtain 

'^=' l,k 

1 < {/«} < q 
- eg J] w,-;,^Lii+,,,,,.. (AJ-'^ • ■ • , A^---'-, Af e) 

1 < {/a} < g 

a=l ? -, 

J, S 

1 < {jm} < q 

E ^/,|Li| (A^^ Aj-'3^ ■ ■ ■ , Aj—^™, A^,™e) , (33) 

J,k 
1 < {/a} < q 

and the next iteration produces 

^^■+i(0 = - E \""' E ^jX^Kx (K~''^ ■ ■ ■ ' ^^,^~^^ K^) 

a=l ? -> 

J, S 

1 < {jm} < q 

E «^^ia (^^^"'^ ■ ■ ■ ' K~''^'^ K^'O ■ m 

J, s 

1 < {Jm} < q 

In this way, Lemma I is seen to be vahd also at order j ' + 1. Consequently, Lemma I is 
valid at arbitrary order. 
Remark I: 

There is another possible parametrization of Eq. (|22|) . It is possible to consider the 
variable C, = 1/C, instead of the variable ^, so that 

^ (35) 



11 



In terms of the new variable, Eq. (122|) takes the form 
1 d _ ~ l^-p-i 



(36a) 



1 ^ ~ (^^ 



fai + 02) :^ + 



1-e^ 



e 






(36b) 



The result of Lemma I does not change under such a reparametrization. 
Remark II: 

In the region < z < 1, the variable ^ is purely imaginary. It is then possible to 
introduce a new variable, y = {1 — 0/(1 + C)i which parameterizes the complex unit 
circle, so that y = exp{i6). The trigonometric parametrization can be derived by putting 
z = tan'^(^/2)/[l + tan''(^/2)], and, for g = 2, it coincides with the parametrization of 
Ref . [32j . In this region, the multiple polylogarithms can be split into real and imaginary 
parts as in the case of the classical polylogarithms [57]. A few particular values of multiple 
polylogarithms, for g = 2, 6 and z = 1/4 (6* = 7r/3), were evaluated in Refs. [32|34II48|I19] . 

2.1.3 One upper parameter is a rational number 

Let us analyze a Gauss hypergeometric function where one of the upper parameters is a 
rational number, 



uj(z) = 2-^1 ( - + «i^, 0-2^; 1 + ce;z 
\(1 



(37) 



Using the algebraic relation between Gauss hypergeometric functions of arguments z and 

1 - 1/z, 



2F1 



a, b 



and putting 



\ 

z 

) 


1 r(c)r(c -a-h) I' a,l + a-c 
~ 2«r(c-a)r(c-6) ^ ^\l + a + b-c 


1\ 

1-- 








, ^a-cn ^^c-a-fcr(c)r(a + 6-c) fc-a,l-a 
^' ^^ '> r(a)r(6) '^'[l + c-a-b 


1- 


1 

z 




P 

a = a2e , b = - + aie , c = 1 + ce 


■> 







(St 
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we obtain 



l+ce 



r(i + c£)r 1 



2 + (c-ai -a2)e 



2Fl 



(a2 — c)e, a2e 



+ 



n _ 2U-p/g+(c-ai-a2)e^(a2-c)e r(l + Ce)! f | + (Oi + 02 



1 
1 - - 

Z 



^c- a2)e 

(c- 02)5, -026 



r(l + a2£)r(£ + ai£ 



X 2;-— 2-ri I -, iD 

rfz VI- f +lc-ai -a2j£ 



1 
1-- 

z 



(39) 



According to Lemma I, the Gauss hypergeonietric functions on the r.h.s. of Eq. (!39|) are 
expressible in terms of multiple polylogarithms with arguments being powers of g-roots 
of unity and the variable r defined as 



r = el.^i-i 



[l-z)-^ 



In terms of this variable, we have 



11(1 - A^, 



T 



:i-z 



.he 



Me 



(40) 



(41) 



The e expansion of the first factor only produces powers of the logarithms ln(l — A^r), 
whereas the e expansion of the second factor generates powers of ln(r). In this way, we 
obtain 
Lemma II: 

The all- order e expansion of the Gauss hypergeometric function 2F1 i-+aie,a2e;l + ce; z) 

is expressible in terms of multiple polylogarithms times powers oflnr, where the variable 

T is defined by Eq. Ij^U^, whereby the arguments of the multiple polylogarithms are powers 

of q -roots of unity times t. 

Remark III 

For «! = 0, all powers of logarithms are factorized, so that the result is expressible just 

in terms of multiple polylogarithms. 



2.1.4 One upper and one lower parameter are equal rational numbers (zero- 
balance case) 

In a similar manner, we can study the so-called zero-balance case [36]. Using the trans- 
formation z -^ —z/{l — z), 



iFi 



a, b 



1 T^ (a,c-b 

2-ri 



[l-zY 



\-z 



(42) 
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and putting 



a = aiE, 



1 \-a2e, c=l \- ce 



we obtain 



1 



{i-zY 



-aie Tp 
2-^1 



aiE, [c—a2)s 



(43) 



According to Lemma I, the Gauss hypergeometric function on the r.h.s. of Eq. (143|) is 
expressible in terms of multiple polylogarithms with arguments being powers of g-roots 
of unity and the variable tj defined as 



zi 



in agreement with Ref. 



1 — 2 

In terms of this variable, we have 



(44) 



i_,fe^uU{i-Kv) 



be 



In this way, we obtain 
Lemma III: 

The all-order € expansion of the Gauss hypergeometric function 2F1 ( 1 — - + aie, a2e; 1 — - + ce; z 

is expressible in terms of multiple polylogarithms with arguments being powers of q-roots 
of unity times the variable tj defined by Eq. ( f^^P - 



2.1.5 All three parameters are equal and non-integer 

In order to derive the e expansion for a Gauss hypergeometric function with three rational 
numbers, the following relation can be applied: 



2-^1 



'1-Z 



c—a—b 



'Fi 



c — a,c — b 



(45) 



Putting 



we obtain 



p p p 

a = l \-ae , 6 = 1 \-be , c=l \-ce , 

q q q 



,l-2 + a£,l-2 + te 



(1 — 2:Y-p/q-{c-a-b)£ 



iFi 



(c — a)e, (c — b)e 

l-2 + C£ 



(46) 



The r.h.s. of Eq. (14611 is expressible just in terms of multiple polylogarithms. In this 
way, we obtain 
Lemma IV: 

The all-order e expansion of the Gauss hypergeometric function 2-F1 ( l — --\-aie, l — --\-a2e] l — --\-ce; z 

is expressible in terms of multiple polylogarithms with arguments being powers of q-roots 
of unity times the variable ^ defined by Eq. ( f7g|) . 
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2.2 The Laurent e expansion of Gauss hyper geometric functions 
with rational values of parameters around 2^ = 1 

Let us construct the iterative solution of the differential equation ([5]) in the neighbourhood 
of the point z = 1. In accordance with the standard procedure [T21IT9] . we introduce a 
new variable, Z = 1 — 2;, so that the the differential equation becomes 



^(l-^)^-[c-(a+6+l)(l-Z)]^-aM^) = 
and one of its solution is uj{Z) = 2-f'i(a; b; 1 + a + b — c; Z) . Setting 



(47) 



a = aiE 



a2e 



1--+C€ 



we rewrite Eq. f l47j) as 



— fz— 

dZ\ dZ 



P 



+ {ai + a2-c)e ] uj{Z) 



^dZ^""'' 



){ 



Z— + a2£ ) u;{Z) . (48) 



Since the difference 1 —p/q can be written symbolically as r/q, where r < q — 1, Eq. fj48l) 
is equivalent to Eq. flT6l) with appropriate changes of variable and parameters, 



{z, c, p) < > (Z, ai + a2-c,q-p) 



(49) 



so that we can use the results of Section 12.1.21 with appropriate change of notations. 

In particular, the solutions of the differential equations for the functions Pi{Q) and 
uJi{Q) can be written as 



^k{Q) = q 



Pk{Q) 






Pk{t) 



(50a) 



(01 + 02) 



f- 



ti 



ai + a2 — c 



t 



Pk~i{t) - qaia2 



Q 



tP- 



ti 



^k~2{t) 



(50b) 



where the new variable Q is defined as 



Q 



z 



1 
1 



(51) 



and ^ is defined by Eq. (TT8l) . 

There is an explicit relation between the solution of the differential equation ([5]) in the 
neighbourhoods of the points z = and z = 1 [19J, which we write in the following form 



iFi 



a, b 



T{c)Tic-a-b) f a,b 
T{c-a)T{c-bf ^ \a+b-c+l 


1-z 
J 


+ ^^ "^ " r(a)r(6) '^'[c-a-b+1 
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1-z 



(52) 



For 



we have 



P 

aiE , b = a2e , c = 1 \- ce , 






r (l - l + c^) r (l - | + (c - ai - a2)e) 



2-ri I p 



'i_zyi~''z~~^'^"~"'~"'^" 



r(l + ai£)r(l + a2e) 






(53) 



The all-order e expansion for the hyper geometric functions entering the r.h.s. of this 
relation is constructed in Section 12.1.21 the l.h.s. is done in this section. Consequently, 
both sides of relation fl53l) are expressible in terms of multiple polylogarithms depending 
on powers of g-roots of unity and the arguments ^ (r.h.s.) and | (l.h.s.). The equality 
of the l.h.s. and r.h.s. of Eq. (!53l) in each order of e generates algebraic relations between 
multiple polylogarithms. 

A similar approach was applied in Refs. [57150] to the connection problem of the formal 
Knizhnik-Zamolodchikov equation in order to derive linear relations between special values 
of multiple polylogarithms. 

In the region < z < 1, the variable ^ is purely imaginary. It can then be rewritten in 
terms of the new variable y = (1— .^)/(l+^), where y = exp{i6). In such a parametrization, 
relation (jSSD is equivalent to algebraic relations between colour zeta values. The algebraic 
relations between particular values of multiple polylogarithms of lower depth and weight 
and particular values of q, specifying the root of unity, were analysed for g = 2 and 
z = 1/4 in Refs. [32]l3lll9] and for g = 6 and z = 1/4 in Ref. [H]. 

3 Multiple (inverse) rational sums 

It is well know that there are three different ways to describe hypergeometric functions: 

(i) as an integral of the Euler or Mellin-Barnes type, 

(ii) by a series whose coefficients satisfy certain recurrence relations, 

(iii) as a solution of a system of differential or difference equations (holonomic approach). 

For functions of a single variable, all of these representations are equivalent, but some 
properties of the function may be more evident in one representation than in another. 
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In Section I2.H the third approach, the iterative solution of differential equations, was 
used to construct the all-order e expansion of a Gauss hypergeometric function. Now, we 
wish to analyze the series generated by the e expansion of a generalized hypergeometric 
function with one rational parameter. This was properly analyzed for the zero-balance 
case in Ref. [36] and for g = 2 in Ref. [39] . 

3.1 Gauss hypergeometric function as generating function of 
multiple (inverse) rational sums 

The starting point of our consideration is the Taylor expansion of the F function. The 
proper expression may be extracted from Refs. [T9f36] and reads 



r(/;+i+f+j) V q) ir-j m ^ (,.+/t+i+£+i 



r('l + 2 + ;) ~ , .^m J+'-- , 



{^-i) 



m=l r=l \ ' ' „ 



where \Ef is the psi-function, "^{z) = ^lnr(z), k is an arbitrary non-negative integer, 
A; > 0, and we have used the two auxiliary expressions 

(7 \ m oo ^ 



r))m+l ' 

n ^ 
^(1 + 2 + n) - ^(1 + z) = V . 

i-^k + z 
fe=i 

In particular, for p = 0, we have 

i"^WT^ = '"r(i + i)-E^^5„.W, (55) 



r(l + z) ^-^ m 

where S'a(j) is the harmonic sum defined as Sai^f) = J2k=i 'h- Based on previous experi- 
ence [52]|34|45j . we choose uj{z) = 2Fi{l + aie, l + a2e; 2 — - + ce; z) as the basis function, 
appearing on the r.h.s. of Eq.®. Using representation ([6]) and expression flMl) for each F 
function, we write the e expansion of our basis function as a multiple series, acl 

(l + a,eA + a,e \ 1 A P^ ^^ /^^'^^ (^~') . ,.., 

2^1 9 P+rP ^]=~\^ — +ce\^z^ — - — ^ ^A, 56 



"^The expansions of T functions about rational values of their parameters may be rewritten in terms 
of multiple Z and S sums (for details, see Ref. [36]). 
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oo / N J. 



E -f- (-^^^^(^' - ^) + "'^!^"'''^ (^' - ^0 



where 

A = exp 

.fc=i 

with Afc = a^ + Cg. Here, 

r=l ( r + £ 

denotes the generahzed muhiple harmonic sum, which satisfies 

1 



k ' 



jb.-j] 



b.e] 



sro+i) = sr'o)+- 

The harmonic sum Sa{j) is a special case of S^' (j), for p = 0, 
In particular, the first few coefficients of the e expansion read: 



h r(i+j- 



exp 



oo / \ I, 









A,S^ij-l)-cSt'''\j-l 



AlSfU-l) - A2S2U-I) - 2cA,S,{j-l)St'''^ (j - 1) 



+c' 



[q-p,q] 



sr"'"' u - 1) 



4q~p,q] 



(j - 1) 



0(e' 



From the relation 



2^1(0, b; c; z) = —2Fi{a + 1, 6 + 1; c + 1; z) , 



dz 



c 



it follows that 








\^ i-er 


[h k 


00 
i— n 


aia2 


e', 





-A,^.(j-l) + c'^Si™l(j-l 



(57) 



(5^ 



(59) 



(60) 



(61) 






(62) 
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where ^ and pk{0 ^^^ defined by Eqs. (fTSll and (120|) . respectively. The algorithms for 
the analytical evaluation of Pk{0 were presented in Section I2.1.2I The equality of the 
l.h.s. and r.h.s. of Eq. (|62l) in each order of e allows one to express the combination of 
generalized multiple rational sums in terms of multiple polylogarithms. Using Eq. ( 13T1) . 
we obtain 



oo r f 1 - 2) r(j) 1 



.=1 r(i-|+J 



ii=i 



^ .r(i-|)r(j) 






(j-i) = e E K'^h,i{K'''^KO 



jl,J2 = l 






.=1 r(i+i-f 



- r(i-2)r(j) 



ii=i 






r(i+j- 



Sk-p..] (^- _ 1) 



[q-p,q] 



U - 1) 



-2g2e E K'^'^ iK^) 

°o t(i-At(j) 



le E A^,^^kil,2(A^^-^^eAf)+Li2,l(A^;-^^eAf) 



iij2=i 



r(i+j-f 



(63a) 



(63b) 



(63c) 



(63d) 



(63e) 






-e E ^ri 2Liu.i (A^^-^^Af-^^^Af) - Aj'^-^^^Lii,!,! (A^-^S A^-^^eAf ) 
iij2,n=i 



fs, .r(i-f)r(j) 



(63f) 



E 



r 1+j- 



(j-i) = e E A(^■l+^■^-^■^Hil,l,l(AJ-^^Af-^^eAf) 



ilJ2J3 = l 



(63^ 
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where ^ is defined by Eq. (TTSj) . Symbolically, Eq. (162!) may be written as 

- ^, r f 1 - 2) r(j) 



pr^^r 1-2+, 






(64) 



J,k 

1 < {jm} < q 

Yll=i ^i = l + c+ai + a2 + . . . + ap 



where A; is a multi-index, k = {ki,- ■ ■ ,kr}, v^ are rational numbers and u^ are complex 
numbers (m^ G C). Unfortunately, we cannot treat all multiple inverse rational sums 
using the all-order e expansion of Gauss hypergeometric functions, but just certain linear 
combinations (see the discussion in Ref. [32]). However, for the evaluation of the others 
another technique can be used. 



3.2 Multiple inverse rational sums of arbitrary depth and weight 

There is an important subclass of multiple inverse rational sums, which are defined as 



- ^,r(j)rfi-£' 

;-W = E - / ^ < S,AJ - ^)SaAj - 1) ■ • • SaAj " 1) , (65) 



t^ r 1-^+j 



-'ai,--- ,afe; -;c 



where ai, • • ■ ,ak,c are arbitrary positive integers. The number w = c + l + ai + --- + ak 
is called the weight and d = k the depth of the sums. For the analysis of these sums, the 
generating function approach P^39II42|[M ] can be applied. 



Let us rewrite the multiple sum fl5^ in the form S|?: 



E^l 2:^ ^a;-;c;-(j) 



where a = (ai, . . . , Op) denotes the collective list of indices and rjs-_.c\-{j) is the coefficient 
of z^ . In order to find the differential equation for generating functions of multiple sums, 
it is necessary to find a recurrence relation for the coefficients r]^._.^._{j) with respect to 



the summation index j. Using the explicit form of i]^ 
the coefficients can be written in the form 



[P.9] 
a; — ;c; 



_{j), the recurrence relation for 






where the "remainder" rg_._[j) is given by 

r(i+j- 



(66) 



r(j)r(i-f 



rt:%)=J^{\{ 



.r=l 



J . 



n^-(^'-i) 



(67) 



r=l 
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Multiplying both sides of Eq. (l66l) by z^ , summing over j = 1, 2, 3, . . ., and using the fact 
that any extra power of j corresponds to the derivative z{d/dz) leads to the following 
differential equation for the generating functions E|? .^._(z): 



z J dz z q 



-^) ^;tM) = s,,o + Bi;^{z), (6^ 



where the non-homogeneous term R^. (z) = Xlfli ^"''"!^_ (j) is again expressible in terms 

of sums of the same type, S{f' . ^ . _.„._(z), but with smaller depth, and 6a,b is the Kro- 
necker 6 symbol. The boundary conditions for any of these sums and their derivatives 
are 

{'ly ^'-'^'-'^-'^^^^^ = ^ (j = 0,1,2,...). (69) 

Let us consider the differential equation fl68l) in terms of the variable ^ defined by Eq. flTSl) . 
The notation E^^ (f) will be used for a sum defined by Eq. (I65l) . where the variable z 

a;6;c;— ^ ' ^ ^ 

is rewritten in terms of the variable ^: 

^t.JO ^ ^it..^ (ziO) ^ ^t.Jz) ,., • (70) 



a;b;ci;C2 a;6;ci;C2 ^^ ' ' a;6;ci;C2 



z=z{^) 



In terms of the variable C,, Eq. fl68l) may be split into the sum of two equations. 






From Eq. (ITTaj), it is easy to obtain 



(72) 



or in equivalent form. 



From this representation, we immediately obtain the following lemma, which is a gener- 
alisation of a statement given in Ref. |45j : 
Lemma A 

//, for some integer j , the series S|?_ .(^) is expressible in terms of hyperlogarithms with 

complex coefficients, then this also holds for the sums S|?_ .^j(^) with positive integers i. 

In order to prove Theorem A for multiple inverse rational sums, we prove an auxiliary 
proposition: 
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Proposition A 

For c = 0, the inverse rational sums are expressible in terms of multiple polylogarithms of 
arguments being powers of q-roots of unity and the variable C,, defined by Eq. ( flg|) . with 
complex coefficients Cr^g times a factor ^^,as 

S?;?,.,; _;0;-(^) ^^^ =e E ^/,^^Li, (A^;-^% Af -^% ■ ■ ■ , A>-->, A>e) , 

J, s 

1 < {Jm} < q 

(74) 
where the weights of the l.h.s. and the r.h.s. are equal, i.e. si + ■ ■ ■ + Sr = 1 + Oi + ■ • ■ + a^ . 
Substituting expression (1711) in the r.h.s. of Eq. (172]) . setting c = 1, and performing a 
trivial sphtting of the denominator, we obtain 



S|^!l;-(-) 



z=z{0 



E E V" /^ r^Li, (A^-^^ Af -^^ ■ ■ ■ a>— >, A>t) 

J, s ^-^ ^' 

1 < {jm} < q 

Y, K'' cj;^Lii,, (A^-^^ Af -^3^ . . . ^ A>-->, A>->+s A>+^e) 

J, s 

1 < {jm} < q 



(75) 



In accordance with Lemma A, we have 
Corollary A: 

Fore > 1, the inverse rational sums are expressible in terms of multiple polylogarithms of 
arguments being powers of q-roots of unity and the variable ^, defined by Eq. / flgj) . with 
complex coefficients djg, as 



y[p,g] (^\ 



J, s 
1 < {3m} < q 

(76) 
where the weights of the l.h.s. and the r.h.s. are equal, i.e. si + - ■ -+5^ = l+c+ai + - ■ ■+ap. 
The strategy of the proof of these results is similar to the one adopted in Ref . [39j . We 
reproduce it here for completeness with appropriate modifications. In Eq. (168|) . it is 
necessary to distinguish two cases: (i) i??' (z) = 0,5^,0 = 1, the so-called depth-0 sums, 
and (ii)4P'^l(^)^ 0,(5,^,0 = 0. 
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Let us first consider the multiple inverse rational sums of depth 0, 






In this case, the system of equations (ITT]) has the form 

Q(i - n^^y ^^^fi;c;M) = evL^AO , (rsa) 

We immediately obtain 

^i':^'(O = -E^fLii(A^,0 , (79) 

i=i 

and 

sL^rJ;0;-(e) = -e E ^f Lii (A^^O , (80) 

i=i 

which agrees with Proposition A, and also with Eq. fl63al) . Iteration of the last equation 
produces 

sLf a;40=- E xf-^'^'u,,, {x'.'^wiiO ■ (81) 

k,ji=l 

In accordance with Lemma A, all the next iterations produce results in terms of multiple 
polylogarithms with complex coefficients. For example, 

sL^fJ;2;-(0 = -'? E Af -^'^^^Li^,! (A^^^ A^;0 

k,ji=l 

- E Af -^-^Hii,!,! (AJ-^^ A^,-^^ Af • (82) 

k,ji,J2=l 

Let us analyze the sums of depth 1, 
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The coefficients of the non-homogeneous part are expressible in terms of multiple inverse 
rational sum of depth 0, and Eq. (17T!) takes the form 



For c = 0, the system of equations ( l83l) takes the simplest form, 



^IpM 



JP'Q] 



S!!:'^;0;-(0=e<'^(0, 

Let us ffist consider the case ai = 1. Using Eq. flHOl) . we obtain 

^i;o(0= E ^rLii,! (A^;-^^ Af , 
iij2=i 



(83a) 
(83b) 

(84a) 
(84b) 



^5) 



and the result for Sj? .q._(^) agrees with Proposition A and reproduces the result of 
Eq. (]63bp . For ai > 2, the r.h.s. of Eq. (]84bp is expressible in terms of multiple polyloga- 
rithms with complex coefficients, so that c|t|';_(0 is also expressible in terms of multiple 
polylogarithms with complex coefficients^ Substituting these results in Eq. (I84ap . we 
obtain results in accordance with Corollary A. For c > 1, the desired results follows 
from Lemma A. 

Let us apply mathematical induction. Let us assume that Proposition A is valid for 
multiple inverse rational sums of depth k, 






- .r(j)r(i-f^ 



i=i 



r 1+j- 



s,i<{jm}<q 



z=z{0 

K"~''^K^ , (86) 



where s = {si, ■ ■ ■ , s^}, and Si + ■ ■ ■ Sr = 1 + ai + ■ ■ ■ + a^- Then for c > 1, Corollary A 
also holds for multiple inverse rational sums of depth k, 



y[p,q] (^\ 



E- / A' ^a.(j-l)---^a,(j-l) 



:t^ r 1+j- 



z=z(0 



E djyU, (A^-^^ Af -^^ ■ ■ ■ , A^,^— >, A>e) • (87) 

s,l<{jm.}<q 



^In particular, Eq. ( |63g[ ) may be reproduced easily. 
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For the sum of depth k + 1, the coefficients of the non-homogeneous part may be expressed 
as hnear combinations of sums of depth j (j = 0, ■ ■ ■ , k) with integer coefficients and all 
possible symmetric distributions of the original indices over the terms of the new sums, 
as 



q c?f 



S£?.,a.,,;-;c;-(0=eVSf3.,a.,,;40, (88a) 



J M-p-1 ~ r(i+7)rfi-2) 






r(i+j- 



q 
k 



xV V 1 s..(j-i)---s..(j-i) 

P=0{ii,-,*fe+i)-^^ ^ ^' -J 

where the sum over the indices (ii, • ■ -ik+i) is to be taken over all permutations of the 
list (ai, ■ ■ ■ ,afc+i). If ip+i + ■ ■ -ik+i > 2, the r.h.s. of Eq. (188bp is expressible in terms 
of multiple polylogarithms of weight k with complex coefficients; see Eq. fl87|) . As a 
result of integrating this equation, a^^... a,,^^--iO is also expressible in terms of harmonic 
polylogarithms of weight k + 1 with complex coefficients. 

If ip+i + ■ ■ -ik+i = 1, the r.h.s. of Eq. (]88bp is expressible in terms of multiple poly- 
logarithms of weight k with a common factor ^P; see Eq. ( 186|) . The result of integrating 
this equation is again expressible in terms of multiple polylogarithms of weight k + 1 with 
complex coefficients, as 






i;-(^) = E E /' dt-^djyU, (A^;-^^ • ■ ■ , A>— >, A>0 . (89) 



;i<{im}<gi=l "'^^ " A; 



For c = 0, direct substitution of the previous results into Eq. ( 188aP shows that Proposi- 
tion A is valid for weight k + 1. In this way. Proposition A is proven for all weights. 
Then, for c > 2, Corollary A is also true for multiple inverse rational sums of depth k + 1. 

Applying the differential operator z^ = —-(1 — ^^)^;§ repeatedly / times to the sum 

S|fj' .. a ■ -■c-(-2)) "we can derive results for a similar sum with c < 0. Thus, Theorem A 

is proven for multiple inverse rational sums. 

Remark IV 

For the particular value q = 2 [p = 1), the multiple inverse rational sums (|62|) are reduced 

to multiple inverse binomial sums, which were studied in Refs. [30|32y39j . 



J 

oo 



E^m^-(^-i)---^«.(^-i)- (90) 

j=i J \j) 
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In order to convert the results of Eqs. (163!) and Theorem A to the form presented in 
Refs. [32][39] . respectively, is it necessary to consider the new variable 

1-e 



y 



i + e 



i-e 



In particular, 

SL';1;_(^) = -Lii,i (1, - Lii,i (1, -0 + Lii,i (-1, + Lii,i (-1, -0 



For practical applications, the following relations are useful: 

"y dt 



Lii,i(2;,l/) 



I- z 



ln(l - xt) 



Lii,i(l,2;) = -ln2(l-z), 

Lii,i (-1, z) = ln(2) ln(l - z) - U, (^) + Lis (\ 



(91) 
(92) 

(93) 
(94) 
(95) 



Remark V 

Let us modify the multiple inverse rational sums of Eq. fl65l) by introducing an additional 

parameter d, as 

^ , - ^jT{j)T{d-A 

Equation (!66|) is changed to become 



3 + d-- 



and we have 
1 



-l]z^+'- d-l-P 
z / dz z \ q 



'1^ ^t-:-Ad.^) = 5,,, + Rl:l\d^z) 



■>[v,q\i 



In terms of the variable ^, Eq. ( l98l) is split into two pieces. 



g di 



^-^^,^/'-Sii)=ho+R^'Si) ■ 



(97) 
(98) 

(99a) 
(99b) 



For the analysis of this system, the previous technique can be applied directly after using 
the decomposition 

where d is an integer and qd > p + 1. 
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3.3 Multiple rational sums of arbitrary depth and weight 

Another important class of multiple sums are the so-called multiple rational sums, defined 
as 



^ ai,---,afc; -;c;-l^J 



E 



r(j+ 



W{j + 1)T(1 



Sa,U -^)Sa,ij- I) ■■■Sa.ij- I), (100) 



where Oi, ■ ■ ■ ,afc,c are arbitrary positive integers. The quantum numbers, depth and 
weight, are defined as in the case of multiple inverse rational sums, namely the weight as 
w = c + 1 + ai + ■ ■ ■ + a^ and the depth as d = k. In this case, the proper recurrence 
relation is 

(j + i)^^'4^i;-(j + 1) = {j + -^ f^A.,M) + (j + fj rt:%) , (101) 

where the "remainder" rg._{j) is given by 



r(j+i)r(i+| 



JP.9] 



r J+J 



u) = n 



,r=l 



Sa,.U-l] 



jar 



U^arU-l] 



(102) 



r=l 



[pM 



The differential equations for the generating functions T^, .^._(z) is 



1_A^^_P 
2; J dz q 



.^|TK..W^^.« + 



[p.?] 






where the non-homogeneous term Xl^i ^'^^h--{J) ^^ again expressible in terms of sums of 
the same type, but with smaller depth. Symbolically, we write it as 



i=l s,k 

where dgk are complex coefficients and 



(104) 



Y^a, = k + Y,. 



Introducing the variable r defined in Eq. fHOj) . so that 



2 = 1 - r« 
we obtain 



d _ 1 /l-^^\ d 
dz q\ ti ) dr 



[l-z] 



d p 
dz q 



l(rA+p(^q5) 



l(l-r«) d 

T-— 

q r? dr 



TK._(r)=r-VJ:-](r) 



q dr ' Q' V ''" 



(106a) 
^-p)E^-^tV(-)-(106b) 
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The point z = transforms to the point r = 1, so that the boundary conditions for these 
sums are 

T,-;_;e,;-(l) = . (107) 

From Eq. fllU6ap . it is easy to obtain 



1 (1 - r«) d 
T 

q r^ ar 



c-J 



or in equivalent form, 



1 (1 - r*?) d 

r-- 

q r'J dr 



c-j-l 



a\~\c\—\ ) a; — ;j; — V / 



7-1 



(108) 



■^[p^i] (^\ = -q I ^+ ^ ■v'^P''' 



' ^Y~^^^^^^-^ ^^ - ^^ • ^^°^) 



We wish to point out that the point r = 1 is a regular point of multiple rational sums. 

From representation (I109p . we immediately obtain the following lemma, which is a 
generalization of the statement given in Ref . [15] : 
Lemma B 

If for some integer j , the series T|? ..(r) is expressible in terms of hyperlogarithms with 

complex coefficients, then this also holds for the sums T|?_ ._j_^(r) for positive integers i. 
Let us return to Eq. (]106bp and rewrite it in the form 



lr'-'i-i:^^)=s^,o+j:d,-. 



^Tj-|^(r) + Tj-|^_,(r) 



(110) 



where dg---j is a set of constants. Integrating it by parts, we find 



(111) 



Substituting this expression in the r.h.s. of Eq. O106ap . we obtain 



l(l-r'') d 
T 

q T'i dr 



-b>g] 



q 



-^r^c^-ir) = -^6,^(1 - r^n 



E^^> 



s,J 



tp- 



IpM 



rr,ir) + qr-PJ^ dty^T^.-^i^ 



(112) 



In order to prove Theorem B for rational sums, we first prove the following auxiliary 

proposition: 

Proposition B 

For c = 0, the inverse rational sums are expressible in terms of multiple polylogarithms of 
arguments being powers of q -roots of unity times the variable t, defined by Eq. Ij^U^, with 
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complex coefficients Cr^g cind dr^g times a factor r ^, as 



^ai,-,ap; -;0;-l^J 


z=z{r) 


E {^jy,k + dj,g,,r--)\n'T 


J,s 


< {Jm} < q 





Li.^(A^;- 



-J2 



A>r) - Li, (A^;- 



-J2 



K) 



(113) 
where the weights of the l.h.s. and the r.h.s. are equal, i.e. si + ■ ■ ■ + Sr + k = ai + ■ ■ ■ + ap . 



Substituting expression (11131) in the r.h.s. of Eq. (11091) . setting c = 1, and performing 
a trivial sphtting of the denominator, we obtain 



-[P.9 



T!?f!l;-(-) 



z=z{t) 



(7 , 5, /C 



i=i 



dt 



'-i 



\n^t 



Li,(Ar^^■■■A^'■t)-Li,(A^-^^■■■A^'■) 



1 < {jm} < q 

= - E 

rJ ^ S ^ K 

1 < {jm} < q 



djMr 



Lil,4Af-^^■■■,A>+V)-Lil,4A^- 



■n 



Aro 



(114) 



In accordance with Lemma B, we have 
Corollary B: 

For c > 1, the multiple rational sums are expressible in terms of multiple polylogarithms 
of arguments being powers of q-roots of unity and the variable t, defined by Eq. ( |7^ , with 
complex coefficients dg, as 






z=z{0 



(115) 






%-,fcln'T 



Li, (A^-^^ ■ ■ ■ , A>r) - Li, (A^-^^ ■ ■ ■ , A>) 



(Oil 



1 < {jm} < q 
where the weights of the l.h.s. and the r.h.s. are equal, i.e. Si + - ■ ■ + Sr + k = c+ai + - ■ ■ + ap. 

In order to prove Theorem B for multiple inverse rational sums, it is sufficient to 
proof Lemma B. The strategy of proof of these results is similar to the one adopted in 
Ref. [39]. We reproduce it here for completeness with appropriate modifications. 
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In Eq. f ll06p . it is necessary to distinguish two cases: (i) /??' (z) = 0,6gfl = 1, the 
so-called depth sums, and (ii) Rp (z) ^ 0,6^^ = . 

Let us consider the multiple inverse rational sums of depth 0, 



r(i+l 



In this case, the system of equations (I106P has the form 



Z = 2(t) 



r = r '^ai:::[T 



q r"? dr J ~'~' 
dr 



We immediately obtain 



and 



a!^.lVr) = ^ (1 - r^) 
p 



T^;^;0;-(-) 



Iteration of the last equation produceqj 



p ^ ' 



g-1 



^-';'-;i;-(^) = ! E(l - K') [Lii (A^) - Li^ (A^j 



^7=1 



V 



E(l-^f)[Lii(A^)-Li,(Aa] 



i=i 



where the last term of the sum, for j = g, is identically equal to zero. 
The next iteration yields 

q-\ 
-;-:2;- 



TL^.'i,.Jr) = -^ 5^ (1 - A^;^) [Lii,i (A^;-^^ Af r)-Lii,i (Af-^^ Af )] 



p 



ii.i2=i 

q-\ 



P 



J2 {i-Knu,{x^) 



JlJ2 = l 
q-1 



Lii(Afr)-Lii(Af) 



:il6) 



(117a) 
(117b) 

(118) 

(119) 



(120) 



+^E(l-^r) [Lii,i (r,A^,0-Lii,i (1,A^,0+Li2 {X^r)-U, {X^)] 






[\2\\ 



^We wish to mention that 



g-l 

^hi(l-A^,)=hig. 
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where we have used the identity [9] 

Lim (x) Li„ (y) = Um,n {x, y) + Li„,m (y, x) + him+n {xy) 



:i22) 



In accordance with Lemma B, all the following iterations produce results in terms of 
multiple polylogarithms with complex coefficients. 

Let us now analyze the multiple inverse rational sums of depth 1, 



^ ai; — ;c;—\ ' J 



E 



r(j+ 



V r{i+j)T(i+^, 



.s,^(j-i)^J2 



r(.+?) 



i-i 



r^^r(i+j)r 1+ 



E-- 



i=l 



The coefficients of the non-homogeneous part are expressible in terms of multiple rational 
sums of depth 0, and Eq. (11061) takes the form 



q T'i cLt J 



MM 



For c = 0, the system of equations fll23p read 



dt ^[^'^1 



l-f? 



s'^:^;.,-i(t) . 



T£^;0;-(r) = -T 



'^:^;a,-(r)-gr-^ 



j.p— 1 
dt- S'^''^! 



;«! 



-iW • 



(123a) 
(123b) 

(124) 



Let us first consider the case ai = 1. Using Eqs. (I119p and (I120p . we obtain 



9-1 



Tt;%.^_{r) = -^E [(l-^f) + ^''(l- V')] [Lii (A^,r)-Lii {Xi)]-'Lr~nnT 



P 



i=i 



P 



(125) 



in agreement with Proposition B. For ai > 2, the r.h.s. of Eq. (11241) is expressible in 
terms of multiple polylogarithms with complex coefficients, so that this also holds for 
T|f^'^l(r), in agreement with Corollary B. For c > 1, the desired result follows from 
Lemma B. 

Let us apply mathematical induction and assume that Proposition B is valid for 
multiple inverse rational sums of depth k, 



-^ai,-,afc;-;0;-l'*. 



-E^'- 



r(.+f) 



r(i+j)r 1 



E 

rj ^ S ^ K 

1 < {jrn} < q 



Cj>,. + ^jV,.r-M In'^ r 



5„,(j-l)---5,,(j-l) 



z=z(t) 

Li,(A^-^^■■■,A^^r)-Li4Af-^^ 



> K) 



(126) 
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where s = {si, ■ ■ ■ , Sr} and si + ■ ■ ■ Sr + k = ai + ■ ■ ■ + ak- Then, for c > 1, Corollary B 
also holds for multiple rational sums of depth k, 



^ ai,--- ,ak;-;c;-\^ ) 



1 < {3m] < q 



^^"r(i+j)rfi+?') 



5,,(j-l)---5.,(j-l 



'i J 



dj.Mr 



Li,(A^-^^■■■,A>r)-Li,(Af 
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z=z(t) 

K) 



(127) 



For sums of depth A; + 1, the coefficients of the non- homogeneous part are expressed 
as linear combination of sums of depth j {j = 0, . . . ,k), with complex coefficients and 
all possible distributions of the original indices over the terms of the new sums. Using 
relation flllll) . we obtain 



Ul 



r'^ 



q r? 



E 

S,J 



d 

dr 

-[P.9 



^ ai,--- ,ak+i;-;c;-\' ) 



Cs,-,jrfl\AT) + qr-^d. 



fP-i 
dt '^^^•'' 



l-ti 



Tr.-i(^) 



(128) 



Let us set c = and consider the two cases: (i) J = 1 and (ii) J > 2. For J = 1, the first 
term on the r.h.s. of Eq. (11281) is expressible in terms of multiple polylogarithms. The last 
term on the r.h.s. of Eq. (I128p has the structure of Eq. (I126p . so that, upon integration, 
it will again be expressible in terms of multiple polylogarithms of weight k + 1. Due to 
the common factor r~^, the result agrees with Proposition B. For J > 2, both terms 
on the r.h.s. of Eq. (11281) are expressible in terms of harmonic polylogarithms of weight 
A; + 1, as 



1(1 



d 






EE 

s,J j=i L 



dr 

Cs;-,J 



c-1 



T 



IpM 



ai,--- ,afc + i;-;c: 

dt 



Xr) 



-[p,?] 



t-TV 



r^ri.W + d^^-^J 



dt " 



^-^'^ 



fp-i 



dti 



Tt'l i(ti) 



l_f^-s;-,J 



(129) 



In this way. Proposition B is found to be valid for weight fc + 1. Consequently, Propo- 
sition B is proven for all weights. Therefore, for c > 1, Corollary B is also valid for the 
multiple binomial sums of weight A; + 1. 

Applying the differential operator z4- = —-^—^4- repeatedly / times to the sum 



-[p.?] 



q t1 dr 



^ai ■■ ar- -0(2^)5 "^6 can derive results for a similar sum with c < 0. Thus, Theorem B is 
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proven for multiple rational sums. 

Remark VI 

For the particular value q = 2 {p = 1), the multiple rational sums are reduced to multiple 

binomial sums, which were studied in Refs. 131)13211391 . as 



ai,--- ,ap] — ;c;— \^/ 



E 



rO+^) 






-j-r(|)r(j+i) 

2j\ fZ\3 1 1 



5a,(j-l)---5. (j-1) 



^^(?) 



-5a.(j-l)---5a,(j-l). 



(130) 



In order to convert the results of Eq. (!63l) and Theorem A to the form presented in 
Refs. [3TJI32II39] ■ is it necessary to introduce the new variable x as 



r 



1-X 

1 + X 



4 All-order Laurent expansion of generalized hyper- 
geometric functions with one rational parameter 

In this section, we turn our attention to the proof of Theorem C. It is well known that 
any function pFp_i{a + fh] h + k; z) is expressible in terms of p other functions of the same 
type, as [5^1135] 

p 

— * — * — * ^ — ^ — * — > — * 

Rp+i{a,h, z)pFp_i{a + rh]h + k] z) = } ^ Rkja, b, z)pFp_i{a + el; b + Ek] z) , (131) 

it=i 



where m, /c, e^, and E^ are lists of integers and Rk are polynomials in the parameters a ,b, 
and z. Systematic methods for solving this problem were elaborated in Refs. |52|f53] . For 
generalized hypergeometric functions of Theorem C, let us choose as the basis functions, 
appearing on the r.h.s. of Eq. 013ip . arbitrary p functions from the following set: 



for Eq. fl4bp . r functions of the type 






for Eq. (Hal) , r^ — 1 functions of the type 



p , {l + aisy-'^, {2 + die}'^ 



Q 
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In the framework of the approach developed in Refs. [30]|31ll32f34ll42j . the study of the e 
expansion of basis hypergeometric functions was reduced to the study of multiple (inverse) 
rational sums. It is easy to obtain the following representations: 



r-'^r— 1 



{l + a,£}^, {2 + dieY 
2-£ + fe,{l + e,£}M2 + c,4« 






z 



r^ r—1 



nf=i(i + c^i£)^ r(i 



A 



{l + e,£}M2 + c,£}^ 



inti(i+c.£) 



E 



r(j+ 



Z' 



z u^^^ii+dis) ^ r(i+j)r (i+A J 



K-R-l 



A 



where the superscripts K, L, R, Q indicate the lengths of the parameter lists 



A 



exp 



E 

. fc=i 



k ciP'l], 



wu~' + S,{n - 1)4 + b'SrU) - fSrU - 1) 



pfcob.g]/ 



Sa{n) = Yl^=i l/i" is ^ harmonic sum, and the constants are defined as 



A,^5^af, C.^^cf, Z^.^J^rff, 
tk = Ck + Ek — Ak — Dk, Wk = Ck 



Ek ~- 



E' 



:i32a) 



(132b) 



(133) 



where the summations extend over all possible values of the parameters in Eq. (11321) . In 
this way, for 6 = / = 0, the e expansions of the basis functions in Eq. (I132p are seen to be 
expressible in terms of multiple (inverse) rational sums, which were studied in Section [31 
But all these are expressible in terms of multiple polylogarithms. In this way. Theorem 
C is proven. 



5 Two-loop sunset with two equal masses Af , a third 
mass m, and external momentum q^ = —m? 

The aim of this section is to find a hypergeometric representation for the two-loop sunset- 
type diagrams with the special kinematic configuration considered in Ref. [51] • The 
integral under investigation is 



^122 (">cri, 0-2, m^M^ 




d"A;id"A;2 



[kl-A'PY^[{ki-k2-qY-M'^Y^[kl- 



m 



21a 



(134) 



It is well know that, for arbitrary values of masses and momentum and powers of propaga- 
tors, this integral is expressible in terms of Lauricella functions [25] . The hypergeometric 
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representation of the master integral with three equal masses and an arbitrary value of 
momentum was derived by Tarasov [27] in terms of Gauss hypergeometric functions and 
Appell functions F2. We now demonstrate, that for q^ = —m?, the master integrals 
are expressible in terms of generalized hypergeometric functions with quarter values of 
parameters. 

Using the Mellin-Barnes representation [21], we obtain 



■^122 (« 5 0-1, 0-2, m^M- 



'_-|^ W1+CT2+0+I / j\^^2Nn/2-<Ti-o-2 I 



m 



2\n/2-OL 



-1 /■+ioo 

7-- / du 

2m 



Tia^)Tia2)T{a)T (f ) 

m^y r (cji+u) r {a2+u) r (ai+a2-f +u) r (§ +«) 



n 



T(-u)T{a---u]2Fi 



a — ^ — u, —u 

n 
2 






(135) 



The Gauss hypergeometric function in this expression is of the form 



iF^ 



a, b 
a — b+a 



where, by definition, a is a positive integer. Using the contiguous relations for Gauss 
hypergeometric functions [^Hl^ , it is possible to express this hypergeometric function in 
terms of a linear combination of any two functions with parameters different by integers 
from the original ones. For our analysis, it is sufficient that the proper set of master 
integrals are expressible in terms of integrals with a = 1 [55]. Closing the contour of 
integration at infinity in the right half-plane and summing over the residiues of r{—u) 
and r(l — n/2 — u), we obtain the sum of two hypergeometric series. Using the Kummer 
relation [5611571. 



2F1 



a, b 
1 + a-b 



1 



r(l + a-6)r(l + f) 



r(i + a)r(i + f-6) ' 

where the hypergeometric series on the l.h.s. is defined if a — 6, which is n/2 — 1 in our 
case, is not a negative integer, and substituting it in our series, we obtain a one-fold series 
representation which contains only products of F functions. In contrast to the previously 
studied cases [21], these series contain gamma functions of the type T{a + k/2), where k is 
the index of summation and a is some number. To convert a series of this type into series 
of the generalized- hypergeometric- function type, we split the summation in one over the 
even and one over the odd values of the index, so that 



E/(^)^ 



fc=0 



E 

j=0 



fi2j)hij) + zfi2j + l)h[j 



.2V 



(136) 
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where h and / are arbitrary functions[] Returning to our series, we like to mention that, 
after sphtting the summation into even and odd parts, the term l/F (^-^), which comes 
from the residiues of r{—u), only contributes to the even part. Finally, we obtain 



j[J{l,ai,a2,m^,M^ 



0-1+0-2 ( /\/2\n/2-cri-CT2 (,^2\n/2-l 



(-1)'^i+'^2(M2) 



m 



T{ai)T{a2) 

r(i-f)r(ai)r(a2)r(ai+a2-f) 

1(0-1+0-2) 



6-^5 



ai cri+l £2 CT2 + 1 2cri+2g-2— n 2o-i+2cr2+2— n 
2 ' 2 ' 2 ' 2 ' 4 ' 4 

n (71 +(72 g'i+a-2+1 g'i+g-2+2 (71+0-2+3 
2' 4 ' 4 ' 4 ' 4 

(2-n)/2 



AM^J 



rrl_\ ''-''"'' V (1+ai-f) r (1+^2-1) r (l+ori+a2-n)r (f-1) 

^' M^J r(f)r(ai + a2 + 2-n) 



-| 2+2(71— n 2+2(72— n 4+2(7i— n 4+2(72— ra o'i+(72+l— n (7i+(72+2— n 

'4'4'4'4' 2 ' 2 

T-'^d \ 2+n 6— n (7i+(72+2— n (7i+(72+3— ra (7i+(72+4— n (7i+(72+5— n 

4 ' 4 ' 4 ' 4 ' 4 ' 4 

(4-n)/2 



m 

4M4 



n(4-r2) V^V 



/m^\ ^"-"^/' r(2+ai- f) r (2+a2-f ) T (^1 + ^2 + 2-^) 



r(ai + a2+4-n) 



7-^6 



-| 2(71+4— ra 2(72+4— ra 2(7i+6— ra 2(72+6— n (7i+(72+2— ra (7i+(72+3— ra 

-'^'4'4'4'4' 2 ' 2 

4+ra 8— ra (7i+(72+4— ra (7i+(72+5— ra (7i+(72+6— n (7i+(72+7— n 

4 ' 4 ' 4 ' 4 ' 4 ' 4 



4M4y 



(138) 



Closing the contour of integration at infinity in the left half-plane and summing over 
the residiues of the F functions, we obtain hypergeometric series in terms of the variable 
M'^ jw? . The same results can be derived from the general formula for the analytical 
continuation of the generalized hypergeometric function. 

The first and the last hypergeometric functions in Eq. (11381) may be reduced by the 
Takayama-Zeilberger [521153] algorithm to the basis function 



4-^3 



l+flie, 1+02^, I + cie, I + C2e 



2+de,2-l + bie,2-l + b2e 



(139) 



and its derivatives. The second hypergeometric function in Eq. (11381) belongs to the class 



4-^3 



l + ttiE, 1 + 02^, 1 + 036, I + C£ 

| + rfe,2-i + 6i£,2-| + 62£ 



z 



(140) 



^ In particular, for the hypergeometric function, we have [56158] 



,ar 






,bq 



Qp Qp + l 



/ (U ai + l 
TP I 2 ' 2 ' ' ' 2 ' 2 

2p^2,;+l 1 f,^ fci + 1 h, bg + 1 



nL 



a, 



nj.r&. 



2 ' 2 ' 2 ' ' 2 ' 2 

/ ai+1 ai+2 

P I 2 ' 2 ' ' ' 

2pJ"2(j+l I 3 fc^-i-i fc^+2 

\ 2 ' 2 ' 2 ' ' 



gp + l ap+2 
2 ' 2 

fcg + 1 bg+2 

' 2 ' 2 



4^-9-1^2 



(137) 
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The 8 expansions of the hypergeometric functions appearing in Eq. (I138p are expressible 
in terms of elhptic integrals, which have been found in Ref . [51] , and their generahzations, 
which are beyond our consideration (see the discussion in Ref. [45]). Three master- 
integrals in n = 4 — 2e dimensions, where e is a parameter of dimension regularization 
[2], correspond to cxi = (T2 = 1, o"i = 1, (T2, o"i = o"2 = 2. We obtain from Eq. (I138p 



4-^(1,1,1, m',M') = 


-- - 


-(M^)-(m^)^-I 


2(l+£ 
^{1-8 


) 






p -^5 2' 2' 2 ^ 2 
V ^-^'4' 4 


4MV 


/M^V^' 1 / 1 _! + £ £ e 

1 I? -^'2^2'2'^ 

Uv (i-2e)' 'Vl-|,i + i,! + i 


4MV 


(M\' (l-e) / l,i + f,i + e,e 
\m^J (2-8)il + 28)' 'V2-f4 + l'i + l 


4MV ' 


41(l,l,2,m^M2) = liMr'-'im'r-^^^^}^ 




/l i l + £ l + £ 
„ -^5 2'2'22 
4-^3 f) ^35 


4MV 


fM^V^'l ( 1 1+e ^ 8 

I 2 ) c-^^ I 3 e 1 I e 3 I e 
\Tfl J 8 \2 2'4~'~2'4~'~2 


4MV 


Vm^j (2-e)(l + 2e)' ' V 2 - f , | + 1' ! + 1 


4MV 


jy (1, 2, 2, m\ M') = -{M'r'-%m'y-^T\l + ^)^|^ 




1 /l 3 ^ , £ 3 , £ 
P '2' ^2' 2^2 

V -^ ~ ^5 4' 4 


4MV 


fM^y-' 8 /i,i+e,i+f,i+£ 

Vm^j (i + e){l + 28f 'V|-i4 + l4 + l 


4MV 




/M2\-" (1- 
Vm^J (2-£) 


- 8 

;3 


1 


^^14 + 
2.)^'^^U-| 


|, l+e, 1+e 

1 + £ Z + £ 
' 4 ^ 2' 4 ^ 2 


4MV 



;i41a) 



(141b) 



'141c) 



To cross-check our resuhs, we evaluate the first few coefficients of the expansion of the 
original diagram in the large-mass limit [59] using our program packages [60] and compare 
it with the proper 8 expansion following from the hypergeometric representation (11411) . 
In the latter case, using the relation 



F 1 



{ai} 



p~ii p p ^ 



n-|6, 



l,{l + ai} 
2, {1 + 6,} 



(142) 
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all hypergeometric functions are reduced to one of the two "basic" ones, 

ki + l + aie,k2 + l+a2e,ri + ^ + Ci£:,r2 + | + 026 



^ ^' k3 + 2+de,r3+l + l + bie,u + l + l + b2e 
ki + l + aie,k2 + l+a2e,ks + l+a3e,ri + ^ + Cie 



4-^3 



r2 + l + C2e,r3 + l + ^ + bie,r4 + l + l + b2e 



(143) 



where {ka} are non-negative integers and {r^} are integers. Their e expansions can be 
easily constructed with the help of Eq. fl54ap and read 



4-^3 



ki + l + aie,k2 + l + a2e,ri + ^ + Cie,r2 + h + 026 



k3 + 2 + de, r3+l + | + bie,r4 + l + j + 62^ 
1 r(A;3 + 2)r(r3 + l + |)r(r4 + l + |) 
zT{k, + l)T{k2 + l)T (n + i) r (r2 + |) 

^ Y- ^, rik,+j)T{k2+j)T (ri-|+j) r (r2-i+j) ^^^^ 
X (d-Si^'^''^ (j-1) - arS|^-^l (j-1) - a^St'^ (j-T 



m=l 



-e 



m 



4-^3 



^^^g[4.^3,4] (^._i) ^ 6-s|^'-*+^'^i (j-1) - crsl^'--^'^] (j-1) - c™sP^'-i'2i (j-i; 



A;i + 1 + ai£:, k2 + l+a2e,k3 + l + a^e, ri + ^ +ci£: 



r2 + | + C2£:,r3 + l + | + 6i£:,r4 + l + | + 62£ 
i r(r2 + |)r(r3 + 2-l)r(r4 + 2-|) 

z r(fci+i)r(A;2+i)r(A;3+i)r (n+l) 



X 



E 



r(fci+i)r(A2+i)r(As+j)r (ri-i+j) 



exp 



,=1 r(j)r (r2-i+i) r (r3+|+j) r (r4+i+j) 

X f-aTSt'^ (j-1) - a™Si^-il (j-1) - a-S|t-i] (^_i) 



m=l 



(-^)' 



m 



+^^g[4.^3,4] (^_i) ^ 6-s|:^^*+^'^l (j-1) - crS^r^'^l (j-1) + c-S^I^'^-i'^l (j-1) 

(144) 

where Sp ' (j) is defined by Eq. ( !58l) . Both kinds of expansions, the large-mass and hy- 
pergeometric one, yield the same results, which also agree with Eqs. (16)-(18) of Ref. [51] . 

6 Discussion and Conclusion 

The proof of Theorems I includes two steps: (i) the algebraic reduction of Gauss hy- 
pergeometric functions of the type in Theorem I to basic functions and (ii) the iterative 
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algorithms for calculating the analytical coefficients of the e expansions of the latter. Step 
(i) is well known [I^PSJ . and, in step (ii), the algorithm is constructed for rational values 
of the parameters [see Eqs. (!28ll . (J39ll . (l43l) and ( l46l) ]. This allows us to calculate the 
coefficients directly, without reference to multiple sums. It is interesting to note that the 
Laurent expansions of the Gauss hypergeometric functions with one rational upper pa- 
rameter are expressible in terms of multiple polylogarithms times powers of a logarithm, 
as shown in Section 12.1.31 We presented in Eq. fl53l) an algebraic relation between Gauss 
hypergeometric functions which allows us to find linear relations between special values 
of multiple polylogarithms. 

We constructed an iterative solution for multiple (inverse) rational sums defined by 
Eqs. (!65l) and fllOOp . It was shown that, by appropriate change of variables, defined by 
Eqs. (fT8|) and (HOj) . the multiple (inverse) rational sums may be converted into multiple 
polylogarithms (see Theorem A and Theorem B). Symbolically, this may be expressed 
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as 



^ .r(j)r(i-2 

E ^' / ^ ^^^^3 - i)5.,(j - 1) ■ • ■ 5,,(j - 1; 



.=1 r(^i+J-f 

= e E 

J, s 

1 < {3m} < q 

Y7k=i Sk = l+aiH http 



CjMs{K'''^K 



2-J3 . 



z=z{0 



E - / ^ < s^.U - i)5a.(j - 1) ■ ■ ■ s^,{j - 1) 



i=i 






2 = 2(0 



/^^ T 1^ ( \^^~^^ \i2-J3 



q 5 5 



J, S 

1 < {jm} < q 
X]fc=i Sfc = 1 + c+aiH h( 



.=1 r(j + i)r(i+ 



Sa.ij - l)Sa^{j - 1) ■ ■ • S'a,(j - V 



z=z(t) 



E 



r(j+ 



(145) 



^>-i->^^>^) (c>l),(146) 



Lis (A^^ 



-^^ ■■■ A>r)-Li,(A^-^^ 



q 1 5'>q 



g ' ' 'g 



A^) 



(147) 



r^ r(j + i)r 1+ 



Sa, (j - 1)5^2 (j - 1) ■ ■ ■ Sa^ (j - i; 



U ^ S ^ rZ 



T 



Li,(A^-^^■■■,A>r)-Li,(Af 



-J2 



2=2;{t) 



(C>1) 



:i48) 



where {cj^, cjg, dp^g, dp^} G C are numerical coefficients, the weight of the l.h.s. equals 
the weight of the r.h.s., and 



i=l 



is a harmonic series. 

Unfortunately, one of the unsolved problems is the completeness of the representation 
in Eqs. (11461) and (11481) . In other words, is it possible to express all multiple polylogarithms 
in terms of multiple (inverse) harmonic sums? If not, what kind of sums must be added 
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to obtain a complete basis? Another problem beyond our present considerations is to find 
the algebraic relations between the sums. 

Using the results of Theorem A and Theorem B, we proved Theorem C about the 
all-order e expansion of a special class of hypergeometric functions. The proof includes 
two steps: (i) the algebraic reduction of generalized hypergeometric functions of the type 
specified in Theorem C to basic functions and (ii) the algorithms for calculating the 
analytical coefficients of the e expansions of the basic hypergeometric functions. The 
implementation of step (i), the reduction algorithm, is based on general considerations 
made in Refs. [52|53] . In step (ii), the algorithm is based on the series representation of the 
basis hypergeometric functions defined by Eq. f ll32p . The coefficients of the e expansions 
are expressible in terms of multiple (inverse) rational sums, to which Theorem A and 
Theorem B apply. 

Finally, we demonstrated, in Section [5l that Feynman diagrams produce hypergeo- 
metric functions with quarter values of parameters. 
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A Hyperlogarithms and multiple polylogarithms 

For completeness, we present a definition of multiple polylogarithms and some rela- 
tions that are useful for our considerations. We are guided by the analyses presented 
in Refs. piTnll6T]l62] . 

The starting point of our consideration is the integral 



lyO'k-, Ofc-i ■ ■ ■ , cti; 2) 



dtk r" dtk^i _ _ n dh 

tk — O'k Jo tk-i — ak-i Jo ti — ai 

dt 

/(afc_i--- ,ai;t) , (149) 



t — CLk 



where we assume that all a^ 7^ 0. In early considerations by Kummer, Poincare, and 
Lappo-Danilevky [7118] . this integral was called hyperlogarithm. It was treated as an 
analytical function in the single variable z, the upper limit of integration. Goncharov [9] 
analyzed it as multivalued analytical function of ai, ■ ■ ■ , a^, z. One of the properties of 
hyperlogarithms is the scaling invariance, 

/(ai,---,afc;^) = j(^,---,^;l) . (150) 
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A special case of this integral is the following one: 

^mi^,mf^_l,--- ,mi \^ky ' ' ' ; '^li ^J 

= J( 0, ■■■ ,0 ,afe, 0,---,0 ,afc_i,---, 0, ■ ■ ■ , ,ai;z) 
rrife-i times mfe_i-i times mi-i times 

(151) 
where all a^ 7^ 0. The last sum can be rewritten as 

,^,"\?ijr w (J1+J2)- U2, 

Jl — -•- jk — -'- 

1 f z \ ^"-^ 1 / z ^ ^'= 






By definition, the multiple polylogarithm is [D] 

Lifci,fe2,-,fcn (^1)^2, ■ ■ ■ ;^n) = /^ fc- ^ ^, (153) 

m„>r?i„_i>--m2>mi>0 i z 

with weight A; = /ci + ^2 + • • • + fcn and depth n. From relations (I15ip . (11521) . and (11531) . 
we have 

r ( ■ \ — (-A\'^-[- (— ^ — 

*-^m„,m„_i,--- ,m2,m,i y^U) ^n—l) ' ' ' 5 2^2, Xi, Zj (^ ij ijl„ii,m2,--- ,m,„ I ; 1 ' ' ' ) 

\Xi X2 Xn 

(154) 
The inverse relation is 

Lifci,fc2,-,fc„ (1/1, ?/2, • • • jUn) = ( — l)"'G'fc„ fc„_ A:2,fei , , " " " , ] 1 

(155) 
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In particular, we have 



oo 



^ ' ^ Jq t-a ^ja^ ^ a' ^a 

Jo ik Jo 'tk-i Jo h-a ^^j'^a^ \a 



Gi^i{a2,ai;z) = / / 

Jo ti — 02 Jq 



dU 



t2 - fll 

- ) = Li,,, {^,-] . (156) 

7l \ (In / I 7l H- Vo I \ " ~ 

ii=ii2=i 



^^JiV«i/ (J1+J2) va2/ ^'^Vai'«2 



The multiple polylogarithms form two Hopf algebras, the so-called shuffle and stuffle ones. 
The first one is related to the integral representation, the second one to the series. 
Integral (11491) is an iterated Chen integral [63j w.r.t. the differential one-forms 

dy dy 

UJQ = , LL!a = , (157) 

y y-a 



where a is any number, so that 



rrifc — 1 ™fc-i— 1 mi 



mfc,mfc_i,--,miiafc,- ■■ ,«i;^j - / ^0 ^afc^O ^^fc-i " ' ' ^0 ^ai ■ [i-Ofi) 

Jo 
A special consideration is necessary when the last few arguments ak-j,ak-j-i, ■ ■ ■ , a^ in 
the integral I{ai,ak;z) in Eq. (I149P are equal to zero, which is the so-called trailing- 
zero case. It is possible to factorise such a kind of contribution into a product of a 
power of a logarithm and an integral of the type described in Eq. (11491) . An appropriate 
procedure was described for generalized polylogarithms of the square root of unity, the 
Remiddi-Vermaseren functions, in Ref. [^ and extended on the case of hyperlogarithms 
in Ref. [62]. These may be written in the form 

I{A,0;z) = Y,^c^j\nPzh,,k,,-,k„^,{B;z) , (159) 

p,A 

where the coefflcients c ^ are rational numbers and the components of vector B are shuffle 

products of the components of the original vectors A and 0. In Eq. (11591) . the weight of 
the l.h.s. is equal to the weight of the r.h.s. 

B Iterative solution of first-order differential equa- 
tion 

A system of homogeneous linear differential equations, 

^n(t) = A(t,n(t)), 
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where A is an n x n matrix and u is an n-dimensional vector, can be formally written via 
Picard's method of approximation as 

U2 = Uo+ / A{t,ui{t))dt , 

JUQ 

Un = Uo+ A{t,Un^l{t))dt, (160) 

where Uq = u(to) is an initial condition. It can be proven PIM] that, in the region where 
this integral exists, the following properties are satisfied: (i) as n increases indefinitely, 
the sequence of functions Un tends to a limit which is a continuous function of t; (ii) this 
limiting function satisfies the differential equation; (iii) the solution thus defined assumes 
the value uq when t = Iq and is the only continuous solution which does so. 
Let us introduce the set of functions, 

Fp(t) = / rfrA(r)Fp_i(r) , (161) 

J to 

and write the relations O160p as 

u{t) = Fo{t)uo + ■■■ + Rp{t)uo + ■ ■ ■ , (162) 

where Fq is the identical transformation, Fq = L For A{t) = ^ . ^3^, the iterative 
solution coincides with hyperlogarithms of configuration aj. 
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